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Abstract 

The ground states of the one-dimensional Fahcov-Kimball model are investi- 
gated in the small-coupling limit, using nearly degenerate perturbation theory. 
For rational electron and ion densities, respectively equal to |, ^, with p rel- 
atively prime to q and y close enough to ^, we find that in the ground state 
the ion configuration has period q. The situation is analogous to the Peierls 
instability where the usual arguments predict a period-g state that produces 
a gap at the Fermi level and is insulating. However for ^ far enough from 
^, this phase becomes unstable against phase separation. The ground state 
is a mixture of a period-g ionic configuration and an empty (or full) config- 
uration, where both configurations have the same electron density to leading 
order. Combining these new results with those previously obtained for strong 
coupling, it follows that a phase transition occurs in the ground state, as a 
function of the coupling, for ion densities far enough from ^. 
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I. INTRODUCTION 



The theory of the electronic band-structure of solids is one of the oldest theories in 
condensed-matter physics, dating back to 1928 when Bloch proved the existence of electronic 
bands in solidsB. It was soon discovered that most metals could be described with the 
nearly free electron model where the periodic ion potential felt by the electrons is weak. 
Initial emphasis in the development of band theory focused on one-dimensional models 
where both transfer-matrix resultsi and exact solutionJlS could be found. It was only after 
the development of the pseudopotential method! in the 1960's that the success of the nearly 
free electron model was understood. 

In the 1950's, Peierlsi reexamined the perturbation theory for the nearly free electron 
model and found that, in one- dimension, a static distortion would always reduce the energy 
of a solid, because the opening of a gap (at the Fermi level) in the electronic band-structure 
would lower the energy of the occupied electronic states and raise the energy of the unoccu- 
pied states. Such a distortion would produce an insulator from the parent metal and Peierls's 
work led to the conclusion that there can never be a one-dimensional metal. Frohlichi used 
nearly degenerate perturbation theory to show that the decrease in the electronic energy was 
on the order of 6"^ In 6 for a periodic distortion of amplitude 6, while the elastic energy was on 
the order of 5^, so the net effect of the distortion was to reduce the ground-state energy (if 
5 was small enough). Two chemists, Longuet-Higgins and Salemi, independently arrived at 
the same conclusions by examining a general class of Hamiltonians for ring-shaped molecules. 
Note that these conclusions explicitly neglect the quantum fluctuations of the phonons about 
the distorted state, a feature which can lead to a stabilization of the undistorted phase, as 
illustrated in recent rigorous work on the one- dimensional Holstein modehll 

The Peierls distortion is generally studied at half-filling for the electrons, in which case 
the distortion leads to a doubling of the unit cell. The conventional wisdom is that the 
lowest periodic structure that produces a gap at the Fermi level will be the true ground 
state, or, in other words, the Peierls distortion is stable against any higher-order distortions. 
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This result has been recently proven to be true at half fillingllO'B. 

In one dimension, the restriction to half-filling is not necessary, since the same arguments 
used by Peierls in his original discussion, could be extended to arbitrary electron filling, 
with the result that a periodic distortion that produces a gap at the Fermi level will always 
lower the total energy of the system. Furthermore, one can generalize Peierls' arguments to 
include other models, such as models for alloy formation. Historically, these alloy models 
were studied before the models of lattice distortions, where Hume-Rotherylli pointed out 
that certain alloys form only when the electron-atom ratio fell within very narrow ranges. 
In one-dimension, Peierls-like arguments establish the Hume-Rothery rule, since the periodic 
arrangement of the ions will always produce the largest gap at the Fermi level. In the alloy 
picture, it is the periodic ordering of the ions that creates the band structure with a gap at 
the Fermi level, and it is a reasonable assumption to neglect the quantum fiuctuations that 
would be induced if the ions became itinerant since the energetics of interchanging two ionic 
species should be large in a one-dimensional system. 

The simplest model of a binary alloy is the spinless Falicov-Kimball mo deli in which 
itinerant spinless electrons interact with static ions: 



where c] is the creation operator for an electron at site i, Wi is a classical variable that is 
1 if the site is occupied by an A ion and if the site is occupied by a B ion. The hopping 
integral is t and U denotes the difference in on-site energies for an electron on a B site minus 
that of an electron on an A site. The electron density, pe, is the number of electrons per 
site (pe = < c-Ci >), and similarly the ion density pi is the number of A ions per site 

{pi = J^i'Wi) with the number of lattice sites. We are interested in the thermodynamic 
limit, and so we take N (yo, but maintain finite values for both pe and pi. In the alloy 
picture {U > 0), one can envision that the electrons are donated by one of the ionic species 
(say the A ion), in which case a study of the neutral case (where pi = pe) becomes most 
relevant. 
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This model was first proposed by Falicov and KimbalH to discuss semiconductor- 
insulator transitions in certain rare-earth oxides and borides. At zero temperature, the 
electrons are localized on the sites of the metallic ions. As the temperature is increased, itin- 
erant electrons and localized holes are simultaneously produced by excitation. The Coulomb 
interaction is screened out, except for the on-site hole-hole repulsion, which is very strong 
(hard-core repulsion) and the on-site electron-hole attraction, which was assumed to be re- 
sponsible for the transition. In this case, U is positive, and the system is neutral pe = Pi- 
Later, the same model was employed to study ordering of rare-earth ions in mixed-valence 
systemJlil. In this application, the itinerant electrons are ci-electrons, and the static particles 
are localized /-electrons. The interaction is then repulsive {U < 0), and conservation of the 
total number of electrons leads to Pe + Pi = constant, with the constant usually taken to 
be one. 

More recentljJli, work on the Falicov-Kimball model has focused on another aspect, that 
of periodic crystal formation. In this language the A sites are ions, and the B sites are empty 
sites. The question is whether or not the interaction of the ions with the electrons, coupled 
with the fact that the electrons satisfy the Pauli exclusion principle, is sufficient to cause a 
periodic arrangement of the ions to be the ground-state configuration. This is the language 
that we adopt in this contribution. 

The Falicov-Kimball model is probably the simplest many-body physics problem. The 
many-body-physics aspect of the problem enters when one considers an annealed average 
over all possible ion configurations {wi}. At zero temperature, one searches for the ion 
configuration that minimizes the energy of the electrons subject to the relevant constraints 
that are applied to the system. 

The Hamiltonian in Eq. exhibits two different symmetries^ an ion-occupied-empty- 
site symmetry and an electron-hole symmetry. The first symmetry relates the ground- 
state energy (per site) for the configuration of ions {wi} to the energy of the conjugate 
configuration {w*} := {1 — Wi} 
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Eg,{U,pe,{w:}) = Egs{-U,PeAw,})-Upe , (2) 

while the second symmetry employs the unitary transformation Cj — > (— that changes 
electrons to holes, yielding 

Egs{U,l - pe,{Wi}) = Egs{-U,pe,{Wi}) -Upi . (3) 

These two symmetries allow restriction to the region Pe < | and Pi < \, without a loss in 
generality. 

The Falicov-Kimball model has been actively studied in recent years, ever since Kennedy 
and LiebEl and Brandt and SchmidtS independently proved that the period-two phase is 
the ground state for all U when the electron and the ion densities are both equal to |. Most 
emphasis has concentrated on the one-dimensional model, where numerical studieB indi- 
cated that the system phase separated into the segregated phase (where all the ions cluster 
on one side of the lattice) for large enough interaction strength if pe Pi or pe + pi ^ 1- In 
the other cases, where pe = Pi (the neutral case) and U — oo, or pe + = 1 (the mixed- 
valence case) and U — * — oo, it was conjectured that the most homogeneous phase was the 
ground state. These two conjectures have already been proven to be truellUHl. Another con- 
jecture, based upon the many-body version of Rayleigh-Schroedinger perturbation theory, 
stated that in the small U limit the ground-state configuration will be the configuration that 
produces the largest gap at the Fermi level, and this state was shown to have the smallest 
periodicity that could produce a gap at the Fermi level (consistent with the Peierls picture). 
Recent analytical^ and numericalil work on the neutral case has shown, however, that at 
low electron density, there is a tendency for molecule formation, rather than a homogeneous 
distribution of the ions, and a phase-separated configuration of ions may yield a lower energy 
than a pure periodic phase. 

In the spirit of the nearly free electron model, we establish two rigorous results which 
are valid for U sufficiently small: 

First, we show that if the electron density is pe = | with p relatively prime to g, and 
Pi = ^, with ^ < Pi < for some integer p', then the ground-state configuration is a 



phase-separated mixture of period-g phases, and possibly the empty (or full) lattice. 

The second result is a statement about the stability of the pure period-g phase for Pe = | 
(with p relatively prime to q) and pi = ^. For pi G [pc, 1 — Pc] with pc ~ 0.371 [solution of 
Eq. (^)], the ground state has period q and is the most homogeneous configuration; it also 
has the smallest periodicity needed to produce a gap at the Fermi level. On the other hand, 
if the ion density pj = ^ is smaller than | or greater than |, then the ground state is always 
a phase-separated mixture of a phase with pe = ^, Pi = (or = 1) and a period-g phase 
with pe = 2^ p-' = a rational that is closest to pc in a well-defined sense. For pi = ^ < pc 
or Pi = ^ > 1 — Pc the same is true, i. e. the ground state is a phase-separated mixture, 
except for special values of pi [those satisfying Eq. (|2^)] for which the period-g phase is 
stable. 

These results show that the close analogy with the Peierls instability is valid only for 
Pc < Pi < I — Pc- We view the analogy as follows: For U = (and pc = ^, Pi = ^ fixed) 
any ion configuration is a ground state, i. e. the probability to find an ion at a given site 
is uniform and equals pj. This uniform-density state is the "undistorted state", has no gap 
in the electronic spectrum, and is metallic. For U ^ 0, sufficiently small, a particular ion 
configuration is selected which has period g. It corresponds to the Peierls-Frolich "distorted 
state", which has a gap at the Fermi level and is insulating. For pi < pc ot pi > 1 — pc, the 
ground state is (in general) phase separated and is a mixture of a metallic and an insulating 
state. This situation does not have a counterpart in the standard theory of Peierls and 
Frohlich. 

Finally, the above results establish the existence of a phase transition in the ground state 
of the Falicov-Kimball model when U is varied. For densities such that the ground state 
is a phase-separated mixture (for U sufficiently small) there must be a phase transition 
as U increases. Indeed, for U sufficiently large, the ground state is known to be either 
the most-homogeneous phase or the segregated phase (which is a different phase-separated 
state). 

Our presentation is organized as follows: in Section II the perturbation theory is devel- 
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oped showing the t/^ In U behavior of the ground-state energy for small U; in Section III 
the perturbation-theory results are analyzed to show when pure phases are the ground state 
and when the ground state is phase-separated; a discussion follows in Section IV. 



II. PERTURBATION THEORY 

It is most convenient to rewrite the Falicov-Kimball Hamiltonian in a momentum-space 
representation before developing a perturbation-series expansion for the ground-state energy. 
Using the standard Fourier transform 

1 ^ 

"'^-^Ee '''c, , (4) 
(with the lattice spacing set equal to 1) yields 



H = J2Hk) -UW{0)]alak-UY.W(k-k')alak' , (5) 

k k^k' 

for the Hamiltonian of the Falicov-Kimball model in momentum space. The wave vectors 
k and k' are restricted to the first Brillouin zone (— vr < A; < tt) and e{k) := — 2tcosfc 
is the unperturbed band-structure. W{27rn/Q) is the structure factor of the period-Q ion 
configuration {wi} 

1 Q 

W{27Tn/Q) := ^ E e-'^'""^/'^)^;,. , (6) 

Q j=i 

defined for n = 0,1, ...,Q — 1. (It is notationally simpler here to define the fc-vectors with 
k = 27Tn/Q to sometimes lie outside of the first Brillouin zone. Of course, translation by 
—271 will shift these vectors back into the first Brillouin zone.) Note that W{0) = pi by 
definition. 

We begin by performing the many-body version of Rayleigh-Schroedinger perturbation 
theory with the double-summation term in Eq. (^ acting as the perturbation. The analysis 
is straightforward^, requiring a momentum-space integral that can be evaluated analytically, 
yielding 
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2t 

Egs{U,pe, {wi}) = sm(7rpe) - Upepi 

8TTt^^ sin(7rn/g) 



sin^nn/Q) — sin(7rpe 



+ 0{U') , (7) 



sm{Tm/Q) + sm(7rpe) 
for the ground-state energy of configuration {wi}. 

The perturbative expansion in Eq. (0) has a singularity when the electron density is 
rational Pe = | and the ion configuration has a period that is a multiple of q (with the 
exception of those ion configurations whose relevant structure factor vanishes). It was argued 
heuristically in Ref. ^that the configuration with the maximal singularity (i. e., with the 
maximal value of |iy(27rpe)|) will be the ground-state configuration, and this result agreed 
with the numerical work. However, such logic is flawed, because the expansion in Eq. 
is valid for U /t <^ | In | sin(7m/(5) — sin(7rpe)||, which cannot hold when an integral number 
of electronic subbands are filled (i. e., when pe = |). This result was known by Frohlich^, 
and it arises from the fact that there are degeneracies in the unperturbed wavefunction that 
were neglected in the above analysis. 

It is easiest to see the origin of the degeneracies and how to properly treat them by 
examining the perturbation theory of the single-particle energy levels. Wigner-Brillouin 
perturbation theory is used, because it automatically removes the singularities. The ground- 
state energy is found by simply filling up the lowest available single-particle energy levels in 
the system. These energy levels can be expanded in a perturbation series which yields 

to second order in U . The quasiparticle energy E{k, U, {wi}) appears on both sides of Eq. (H) 
because one must self-consistently solve for the energy in a Wigner-Brillouin perturbation- 
theory expansion. The equivalent Rayleigh-Schroedinger expansion would replace E{k) by 
e{k) in the right hand side of which produces a singularity when k = —Trn/Q because 
e{k) = e{-k). 

At this point, textbooks note that the dominant term in the sum over n, in the right- 



hand side of Eq. (^), is the term where k + is closest to 27r — k, i. e. it is the term with 



n 



closest to Q{1 — k/n). If the other terms are neglected, then Eq. (H) reduces to a quadratic 
equation that can be solved exactly. This procedure is sometimes called nearly degenerate 
perturbation theory because it produces the correct secular equation in the degenerate case. 

However, we choose to proceed in a more rigorous manner in the case where the value of 
the interaction is much smaller than the subband width U <^'Kt/Q. In this case, the effect 
of the additional terms can be treated in a perturbative fashion, which gives 

2-nn U"^ 

E{k, U, {wi}) = -t[cos k + cos{k - -^)] + —fn{k) 



±^{t[cosk-cos{k--^)]-—f4k)y + —\W{-^W + 0{U') , (9) 

with 

I Q-l ||y(2zn21)|2 

fnik) :=-4 E ^Qg^_^^g'^^_2^^ , (10) 

for 7r(n — < k < 7r{n + ^)/Q. The minus sign is for the subband energy with k ^ nn/Q 
from below, and the plus sign is for k 7m/Q from above. This form for the quasiparticle 
energies is exact for all TJ when Q = 2, but is perturbative for all higher periods. 

The ground-state energy is found by summing up all of the quasiparticle energies with 
\k\ < kp = vrpe {kp is the Fermi wavevector). Since the quasiparticle energies reproduce the 
noninteracting result when Tf = 0, the zeroth and first-order terms are correctly produced 
by this summation. We want to concentrate on the higher-order terms. The solution for 
the quasiparticle energies reveals that a generic period-Q configuration will break into Q 
subbands. The band gaps are equal to 2U\W{^^)\ and are symmetrically displaced to 
lowest order; the order f/^ correction leads to asymmetries in the subband-structure. If the 
Fermi energy lies within a subband, then it is easy to show that for TJ '^irt/Q the shift in 
the ground-state energy is of order U^/t, because the square-root in Eq. (|^) can always be 
expanded in a convergent power series in TJ. However, no such perturbation-series expansion 
can be made if the Fermi energy lies within one of the band gaps. In this case the ground- 
state energy actually has a TJ^lnTJ dependence^ which is always larger than any order f/^ 
dependence for small enough TJ. 
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We illustrate the origin of the C/^ In U terms in the expansion for the ground-state energy 
for rational electron densities Pe = | with p relatively prime to q. We consider any ion 
configuration with a period Q that is a multiple of q. This guarantees that there will be a 
band gap at the Fermi momentum kp — T^pe- The ground-state energy is 

Eg,{U, pe, {Wi}) = E E{K U, {Wi}) = - E{k, U, {Wi})dk . (11) 

Since the band gaps are symmetric to lowest order, the effects of the lower filled subbands 
cancel, and the In U contribution arises entirely from filling the uppermost subband. 
Therefore, the U'^ In U contribution comes from the integral 



1,^-1 J{t[cosk-ms{k-2T:p,)]-^U{k)Y + ^\W{2Tipe)\Hk . (12) 

TT Jn{pe-^) V It 



Use of the identity cos A; — cos(A; — 27rpe) = — 2 sinvrpe sin(A; — vrpe) and shifting the integration 
range k — > —k + Trpe yields 



1 rf^Q 

7 = / W[2isin7rpeSinA; fpiirp^ - k)]'^ + U^\W{2Trpe)\^dk . (13) 

TT >/0 V t 

The C/^ In U behavior originates from the region near the origin and fp{npe — k) does not de- 
pend strongly upon k in this region, so we can approximate the integral by replacing sin k 
k and fp{'KPe — k) ^ fp{T:pe). The substitution k — >• {U\W{2TTpe)\ sinhx -|- ^/p)/2tsin7rpe 
yields an integrable form for I which contains a constant term and a U^lnU term. The 
small-C/ expansion for the ground-state energy then becomes 

E,s{U, Pe, {m}) - -- Sin TTPe - U p,p, + ^ M^I^^fl^ [/2 ^ ^(^2) ^ (^4) 

TT Ant sm Trpe 

which contains no fp dependence. The above form is only vahd for U <^ nt/Q. This 
perturbative expansion shows that the ground state will be found by determining the periodic 
configuration {wi} that maximizes the square of the structure factor |iy(27rpe)P at twice 
the Fermi momentum. Furthermore, it eliminates all configurations with periods Q that are 
not multiples of q, since those states only have a f/^ correction to their ground-state energy 
because the Fermi level does not lie within a subband gap. 
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III. PHASE-SEPARATION ANALYSIS 



We are interested in finding the ground state of tlie Falicov-Kimball model as a function 



of the electron and ion densities. The perturbative expansion in Eq. (14) depends on pe 
in the zeroth-order term, which is a convex function of the electron density. Therefore, for 
U = 0, phase separation can only occur between two different ion configurations that have 
the same electron density pe as the pure phase. 

Let us examine the effect of the first-order term. To order U, Eq. (|T^ is a concave 
function of {pe, pi) and thus the ground state will be a mixture of two phases with densities 
(Pe> Pi) and {p'^, p'l). We set pe = ap'^ + (1 - a)p'^, pi = ap\ + (1 - a)p-', and p'^ = pe + Spe- 
Since — sin(7rpe) is convex, then for U = 0, we have 6pe = and the probability that a given 
site is occupied by an ion is pi. Hence, 6pe tends to zero as f/ ^ 0. Furthermore, one can 
check that the minimum of 

aE(p;,p:) + (i-a)E(p:,pn , as) 

(at first order) is attained for 

5pe = .^f , , p:' = , (16) 
2t sm(7rpe) 

and the decrease in the ground-state energy is of the order U"^. This is negligible in compar- 
ison to the In U term so one can assume Spe = at this order. 

It is the coefficient of the U'^ In U term that determines which ion configuration yields the 
lowest energy. Since the electron density is fixed in all candidate ground-state configurations, 
the criterion for selecting the ground-state configuration is to maximize the square of the 
structure factor \W{2'7ipe, {wi})]"^, including the possibility that phase-separated mixtures 
may be needed in the maximization. 

The construction of the maximum square structure factor is a straightforward exercise 
for each phase {wi}. Consider a rational electron density Pe = | with p relatively prime to q 
and a rational ion density pi = — with pi relatively prime to g^. Then the maximum of IWI"^ 
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is achieved with the following period-Q ion configurationllZI (with Q = lcm{g, g^} =: sq): 
Define the q numbers rj by 

{prj) := j modq , j = 0, 1, g - 1 , (17) 

and set Wi = 1 for 

i = Tj + mq , j = 0, 1, n — 1 , m = 0, 1, s — 1 , n = int[gpi] , 
= Tn + mq , m = any {Qpi — sn) numbers in the set {0, 1, s — 1} . (18) 

Note that the above construction is not necessarily unique when s 7^ 1, but every configura- 
tion constructed in such a fashion will have the same square structure factor (the order U"^ 
corrections to the energy should split any remaining degeneracies). It is easy to verify that 
the above construction does satisfy W{0) = pi and 

I a — k) 

|Vr(27rpe,Pi)P = TyiYl WjWkCO?,[2TXpe ] , 

^ j,k=l ^ 

_ n-qpi + {n- qpjf 

II + {qpi -n-l) cos(27rn/g) - {qpi - n) cos[27r(n + l)/q] 

^2g2 sin^vr/g ' ^ ' 

In the special case where s = 1, so that n = qpi, the above form simplifies to 

1^^./ m9 1 (1 — cos27rpj) 1 sin^ TTOj , , , . 

W{27rp,,p,)' = — ^ ■ 2 / = 0^ = 1 • 20 

2g^ sm TT/g g^ sm vr/g 

Note that \W{27ipe, Pi)\'^ depends on pe only through the denominator g. This fact greatly 
simplifies the analysis below. 

The ion configuration that maximizes the square of the structure factor is identical to 
Lemberger's most-homogeneous configuration^ in the neutral case pi = pe- In the nonneu- 
tral cases, the maximal ion configuration satisfies uniform-distribution properties0 in which 
the configuration is composed of clusters of ions, with only islands of size / and l—l appearing. 
Furthermore, these islands are "most-homogeneously" distributed (the most-homogeneous 
configuration is the special case with islands of size 1). 
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Lemma (local convexity of the squared structure factor). Assume the electron density is 
rational pe = ^ with p relatively prime to q, and assume that the ion density is also rational 
Pi = ^ with Pi relatively prime to qi, and satisfies ^ < pi < for some integer p' . Then 
a mixture of ionic phases with ion densities ^ and will have a larger square structure 
factor than the pure phase with ion density pi. 

Proof: We need to show that the maximal square structure factor in Eq. (|l^) is locally 
convex. To do this we must examine the condition for convexity, by computing 

C={p' + 1- qp^i^^Pe, -)P + {qp^ - p')\W{2np,, - \W{27rp,,p,)\' . (21) 

q q 

If C > 0, then the square structure factor is locally convex, and the lemma will have been 
proven. Substituting Eqs. (|19D and ( pO|) into Eq. (|2lD yields C = [qpi — p' + {qpi —p'Y]/q^ 
which is greater than zero for < qpi —p' < 1, which is a condition that holds by hypothesis. 
Q.E.D. 

Comment. The above lemma shows that the search for a maximal square structure factor 
can be limited to those ion configurations that possess the minimal periodicity q needed to 
produce a gap at the Fermi level. It has not determined the global maximum. That search 
will be completed below. The lemma does allow us to immediately prove our first result 
about phase separation in the Falicov-Kimball model. 

Theorem 1 (minimal-period phase separation). If the electron and ion densities satisfy 
the hypothesis of the lemma, then for U sufficiently small (i. e., Uq ^ 1), the ground-state 
configuration is a mixture of two period-q phases, with densities p'i = ^ and p" = ^ (p' or 
p" can be equal to or q). 

Proof: The perturbative analysis of Section II established that in the limit U —>■ the 
ground-state configuration is determined by maximizing the square of the structure factor 
evaluated at twice the Fermi wavevector. The above lemma shows that such a search can 
be limited to a search over ion configurations with p[ = ^, p' = 0, ...,q. This means that if 
the ion density does not equal j, then it must phase separate into a mixture of states whose 
electron densities are pp = - and whose ion densities are — and ^. Q.E.D. 
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It follows from the Lemma that we only need to search for the ground state among the 
ion configurations with period q (given Pe = | with p relatively prime to q), therefore s = 1 
and the square structure factor is given by Eq. (pO]). 

The function (cos27rpj — 1) is concave for pi G [0, |] U [|, 1] and convex for pi G [|, |]. 
Hence, if pi lies in the interval [0, ^] U [|, 1], the pure-phase cannot be stable against phase 
separation. 

The convex envelope of the function (cos27rpj — 1) is given by 

(cos27rpc — 1)— for < Pi < pc 

Pc 

(cos 2n Pi - 1) for Pc < Pi < 1 - Pc 
(cos27rpc - 1)^^^ for 1 - Pc < Pi < 1 (22) 

Pc 

where pc ~ 0.3710 is the solution to the equation 

27r Pc = tanvrpc . (23) 

Thus if Pi is a rational in the interval [pc, 1 — Pc], the pure phase with pe = |, pi = ^ is 
stable. 

Let us now analyze what happens for densities pi in the interval [0, pc]. The case [1 — Pc, 1] 
is similar. The Lemma states that we must consider only the ion densities in the discrete 
set {pi = Given pe = |, let ^ be the largest rational in the set {^} which is smaller 
than Pc. From the construction of the convex envelope, for any pi < ^ we know that the 
ground-state configuration is a mixture of the empty configuration p^ = and a period-g 
configuration with density p" = ^ or ^y^. 

To decide between the two possible values of p'-, we have to determine whether ^ corre- 
sponds to a pure phase, or a mixture of the empty state and a period-g configuration with 
density 2i±l. Using Eqs. (|I4|) and it follows that the pure phase — is stable if 

Pi ( Pi \^ . Pi + 1 f^.S 

smvr — > smvr , (24) 

q \pi + lj q ' ^ ' 

and unstable if Eq. (^^ is not satisfied. 
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To summarize, given | and ^ the largest rational with denominator q that is smaller 
than if Eq. ( p4|) is satisfied, then for Pi = ^ (resp. 1 — ^) the ground-state configuration 
is periodic, given by Eq. ([17|) , and for all Pi < ^ the ground state is a mixture with p ■ = 
and p'- = ^ (resp. for all pi > 1 — = 1 and p" = 1 — ^). On the other hand, if 

Eq. (1^) is not satisfied, then for all pi < the ground state is a mixture with p^ = and 
Pi — similarly for pj > 1 — ^. 

In Table 1, we give the values of pi for g = 3 to 34, and indicate whether the pure phase 
with Pi = ^ is stable (s) or unstable (u). For example, the state with pe = ^ is unstable for 
any pj < ^ or p^ > ^ and stable for pj = ^, ^, ^, ^. The state with pe = | is unstable 
for Pi < I or Pi > I and stable for pj = |, |, |, |. In these two examples the neutral state 
Pe = Pi is unstable. On the other hand, for pe = ^ < Pc, the neutral state is stable. 

In general, the neutral state pe = Pi is unstable for pi < pc, with an infinite number of 
exceptions [given by Eq. (p4|)] for which the first few electron densities are pe = |, ^, ^, 
rf' 19' 20' 25' if' M' "^^^ state with diatomic molecules pi = 2pe is unstable for pi < pc with 
an infinite number of exceptions pe = |, ^, ^, ^, ^, ^, ... Similarly, in the triatomic case 
pi = 3pe, the exceptional electronic densities for which the pure state is stable are pe = |, 
Tr' J5' M' • • • ^'^y case, it appears that for any e > and for any state with n- molecules 
Pi = npe there is a finite number of exceptions in [\, Pc — e] as shown in Figure 1. 

These results lead us to the theorem: 

Theorem 2. If the electron density is rational, Pe = g with p relatively prime to q, and 
the ion density is pi = ^, then 

a) for ^ G [pc, 1 — Pc]; or ^ = ^ with pi solving Eq. (^^, the ground-state configuration is 

periodic with period q. 

b) for ^ < Pc (or ^ > 1 — p^) and ^ 7^ ^ with pi solving Eq. (^^, the ground-state 

configuration is a mixture of the empty lattice p^ = and the period-q configuration 
with p'l = (resp. p'^ = 1 and p'l = l-^ ). 

c) for all ^ < ^ or ^ > ^, the ground state configuration is a mixture like in b). 
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Comments: (i) The exceptional ion densities can all be found by studying Eq. (|2^). We 
have not been able to determine an explicit formula for these exceptional ion densities, (ii) 
The phase-separated state is not an insulating state, but rather is the mixture of a metallic 
state (the empty lattice) and an insulating state (the period-g phase with p filled subbands). 
(iii) In the neutral case, Pe = Pi = P, for any p G [pc, 1 — Pc] and for the "exceptional" values 
in the intervals [0,pc] or [1 — pc, 1], the ground state is most homogeneous, since the state 
with the maximal structure factor satisfies the uniform-distribution property. It is also the 
configuration obtained by Lemberger's constructionS. For these pure states it is expected 
that the ground state does not have any phase transition when U increases from +0 to +oo, 
since for any rational density p the ground state is known to be the most homogeneous 
state for U sufficiently large. This expectation is also confirmed for intermediate values of 
U {U > 0.1) by exact numerical calculations^. Using the same argument for the "regular" 
values of p in [0,pc] or [1 — PcPc], there will be a phase transition as U varies, (iv) These 
results only hold for U sufficiently small with respect to ^, where pe = |, for the U^lnU 
term to dominate the perturbation expansion. For pe = ^ and Pi = ^ Pe the phase 
separation that may occur for small U rapidly disappears as U is increased from to oo to 
yield either a pure state or the segregated phaseil. For U sufficiently large, it is expected 
that the state is either neutral or the segregated phase. 



IV. CONCLUSION 

The band theory of solids is perhaps the defining theory for condensed-matter physics. It 
has been applied to virtually every interesting material that has been studied. Nevertheless, 
the conventional wisdom of Peierls and Frohlich for optimizing the band-structure for the 
ground-state of one-dimensional crystals is not always correct. They argue, that the ion 
configuration that produces the largest gap at the Fermi level will yield the ground state. 
We find that this argument is true for a nearly free electron model only if the ion density 
is close enough to half-filling. For ion densities away from half-filling, the system will phase 
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separate into a mixture of states that have the same electron density, but have different ion 
densities (pj = and pi close to 0.371 or pj — 1 and pi close to 0.629). It is possible that 
this phase separation can be observed in quasi-one- dimensional metals and insulators. We 
are not aware of any experiments that have seen this phase separation. 

Our rigorous results hold only for U sufficiently small because they are based on 
perturbation-theory arguments that maximize the leading corrections of the ground-state 
energy as a function of U. Since these corrections of order In U will compete with or- 
der C/^ corrections for finite values of U, the phase separation discovered here may rapidly 
disappear as U increases. Numerical evidence indicates that this is true for the densities 
between | and |, but larger values of U are necessary for the densities near or 1. 

Furthermore, since the ground state is known to be either a different phase-separated 
state (nonneutral cases) or the most-homogeneous state (neutral case) for large U, the 
spinless Falicov-Kimball model must have a phase transition as a function of U. In the 
neutral case, when the ground state is not a phase-separated state, but is the Peierls-type 
state that maximizes the band gap at the Fermi level, it is possible that the ground state 
has no phase transitions for < C/ < -|-oo, since the small- C/ ground state is identical to the 
large-C/ ground state. We are unable to prove this conjecture here. 

Our analysis was restricted to the spinless Falicov-Kimball model, but the general ideas 
may also hold for more complicated models such as tertiary alloy problems (where Wi would 
assume three different values) or the static Holstein model (where Wi is continuous) , but the 
determination of the maximal structure factor becomes much more complicated, since one 
must maximize with respect to both the phase and the amplitude, as opposed to maximizing 
only with respect to the phase, as we did here. 
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TABLES 

TABLE L Largest integer pi such that ^ < pc ~ 0.371. The letters s and u denote whether 
Eq. (24) is satisfied (s), implying the pure phase pi = ^ is stable, or is not satisfied (u), implying 
the pure phase pi = ^ is unstable. 
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FIGURES 

FIG. 1. Stable periodic configuration for n-molecules, i. e. states with pi = npe- The values 

n = 1 (solid dot), n = 2 (open square), n = 3 (solid triangle), n = 4 (open dot), n = 5 (solid 
square), n = 6 (open triangle), and n = 7 (x) are all plotted. The phases are stable above pc as 
indicated by the solid lines. The dashed lines are guides to the eye. 
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